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, We investigate the chiral symmetry breaking and color superconductivity in color and electric 

' charge neutral quark matter including strange quarks (2SC+S quark matter) in a Nambu-Jona- 

' Lasinio model using a variational method. The methodology involves an explicit construction of a 

, variational ground state and minimisation of the thermodynamic potential. Simultaneous existence 

of a mass gap and a superconducting gap is seen in a small window of quark chemical potential 
^ ' within the model when charge neutrality conditions are not imposed. Enforcing color and electric 

\ charge neutrality conditions gives rise to the existence of gapless modes in 2SC+S quark matter 

depending upon the magnitude of the gap and the difference of the chemical potentials of the 
condensing quarks. Within the window for gapless modes, it is seen that the strange quark density 
can be as large as 30%-40% of the number density of either of the condensing light quarks. 
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' The structure of vacuum in Quantum Chromodynamics (QCD) is one of the most interesting 

^ i| questions in strong interaction physics . The evidence for quark and gluon condensates in vacuum 

Ci^' is a reflection of its complex nature 2] , whereas chiral symmetry breaking is an essential feature in 

, the description of the low mass hadron properties. Due to the nonperturbative nature of QCD in 

this regime different effective models have been used to understand the nature of chiral symmetry 



^ . breaking 
5— ( ' 

' Recently there has been a lot of interest in strongly interacting matter at high densities. In 

particular, a color superconducting phase for it involving diquark condensates has been speculated 
with a gap of about 100 MeV at zero tempearture. Possibility of diquark condensates along with 
quark antiquark condensates has been considered in Ref. The natural place to look 

for such a phase seems to be in the interior of compact stellar objects like neutron star. However, 
to apply it to the case of neutron stars, the color and electric charge neutrality conditions need to 
be imposed for the bulk quark matter. Such an attempt has been made in Ref.Q] as well as in 
Ref. where the lighter up and down quarks form color superconducting matter while strange 
quark do not participate in pairing. It has been shown, based u pon comparison of free energy that 
a color superconducting phase would be absent in neutron stars liof. Within NJL model in Ref.[ll| 
it has been argued that such conclusions are consistent except for a small window in density range 



^ Talk presented by A. Mishra at International meeting on Nuclear Physics, Bormio, Italy, Jan26-31, 2004 
'Electronic address: mishra@ th.physik.uni-frankfurt.de, 
t Electronic address: hmlQ)prl.ernet.in| 



2 



where superconducting phase is possible. Novel features of two flavor neutral color superconducting 
matter like gap less modes have been first emphasised in Ref. 123] for zero temperaure and has been 
generalised to finite temperatures [2^. There have also been studies to include the possibility of 
mixed phases [f^ of superconducting matter demanding neutral matter on the average. 

We had applied a different approach to study the problem in Ref. 5] . We considered a variational 
approach with an explicit assumption for the ground state having both quark antiquark and diquark 
condensates. This allowed to simultaneous formation of condensates in quark antiquark as well as 
diquark channels 5] . In the following we generalise this approach to include include the conditions 
of color and electric charge neutrality. This leads to condensate functions which depend upon both 
color and flavor. Although, for simplicity, we shall be considering color superconductivity, this can 
be generalised to the three flavor case to include color flavor locking. In fact here we shall also 
consider three flavors but having the u and d quarks taking part in diquark condensation. Although 
it might look rather complicated at the outset, the niccity of the approach is that within the model 
one can solve for the condensate functions explicitly which are flavor and color dependent j^. 

We organize the paper as follows. In the next section we discuss the ansatz state with quark 
antiquark as well as diquark condensates. In section 3 we consider Nambu Jona-Lasinio model 
Hamiltonian and calculate the expectation value with respect to the ansatz state to compute the 
thermodynamic potential. We minimise the thermodynamic potential to calculate all the ansatz 
functions and the resulting mass as well as superconducting gap equations here. In section 4 we 
give the results and discuss them. Finally we summarise and conclude in section 5. 



2. AN ANSATZ FOR THE GROUND STATE 



For the consideration of chiral symmetry breaking, we denote the perturbative vacuum state 
with chiral symmetry as |0). We shall then assume a specific vacuum realignment which breaks 
chiral symmetry because of interaction. 

The ansatz for vacuum restructuring leading to chiral symmetry breaking that we shall be 
considering is described by 

\vac)^UQ\Q), (1) 

where 

(2) 



Uq = exp (^J q"/{k)^{cT ■ k)hi{-k)q°\k)dk - h.c. 



where, g*^*^ ( qp^) is the two componet quark creation (antiquark creation) operator which acts 
on the vacuum with chiral symmetry denoted by |0). In the above, /ii(k) is a real function of |k| 
which describes vacuum realignment for quarks of a given flavor i. We shall take the condensate 
function /i(k) to be the same for u and d quarks and h^ik) as the chiral condensate function for 
the s-quark. This leads to e.g.(uac|?/5-0|wac) ^ / sin2/ii(fc)dk. Clearly, a nontrivial hiik) shall 
break chiral symmetry having nonzero expectation value for (V'V') ■ Sum over three colors and three 
flavors is understood in the exponent oiUq in Eq.(|2l. 
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Having defined the state as in Eq.lQJ for chiral symmetry breaking, we shall next define the 
state involving diquarks. We note that as per BCS result such a state will be dynamically favored 
if there is an attractive interaction between the quarks Such an interaction exists in QCD in 
the qq color antitriplet, Lorentz scalar and isospin singlet channel. With this in mind, we take the 
ansatz state involving diquarks as 



where 



\n) = Ud\vac) = exp{Bl - Bd)\vac), (3) 



ql^ik)^rr{k)ql'^i~k)h,je3ab + (?T(k)r/r (k)eU-k)e.je3a6 dk. (4) 



In the above, i,j are flavor indices, a, 6 are the color indices and r{= ±1/2) is the spin index. 
As noted earlier we shall have u,d («=1,2) quark condensation. We have also introduced here 
(color, flavor dependent) functions /*'^(k) and fl'^{k) respectively for the diquark and diantiquark 
channels. As may be noted the state constructed in Eq.Q is spin singlet and is antisymmetric in 
color and flavor. Clearly, by construction /'"(k) = /-'^(k) with i j and a ^ b. 

Finally, to include the effect of temperature and density we next write down the state at finite 
temperature and density |r2(/3,/i)) taking a thermal Bogoliubov transformation over the state \ 
using thermofield dynamics (TFD) as described in ref.s We then have, 

\n{p,^,)) ^Upjn) =U0,^UdUQ\o). (5) 

where 1^/3, fi is 

Ufs,^ = e^5n/3,M)-/3(/3./.)^ (6) 

with, 

/ [g;(k)t(?_(k,/3,/.)q;(k)t + g;(k)^?+(k,/3,Ai)g;(k)]dk. (7) 

In Eq.0 the ansatz functions ^±(k, /3, /i) will be related to quark and antiquark distributions and 
the underlined operators are the operators in the extended Hilbert space associated with thermal 
doubling in TFD method. In Eq.Q we have suppressed the color and flavor indices on the quarks 
as well as the functions 9(k, (3, /i). Note that we have a proliferation of functions in the ansatz state 
jrj,/?, /i) - the (flavor dependent) chiral condensate function, the (color flavor dependent) quark 
as well as antiquark condensate functions and the (color flavor dependent) thermal functions. All 
these functions are to be obtained by minimising the thermodynamic potential. This will involve 
an assumption about the effective Hamiltonian. We shall carry out this minimisation in the next 
section. 



3. MINIMISATION OF THERMODYNAMIC POTENTIAL AND GAP EQUATIONS 

We shall work here in a Nambu-JonaLasinio model which is based on relativistic fermions inter- 
acting through local current- current couplings assuming that gluonic degrees of freedom can be 
frozen into point like effective interactions between the quarks. The Hamiltonian is given as 
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Here J° = ip^^T°"ip and, is the current quark mass which we shall take to be nonzero only for 
the case of strange quarks (j=3). Of the two superscripts on the quark operators the first index 'i' 
refers to the flavor index and the second index, a refers to the color index. The point interaction 
produces short distance singularities and to regulate the integrals we shall restrict the phase space 
inside the sphere |p| < A - the ultraviolet cutoff in the NJL model. 

We next write down the energy density which is the expectation value of the Hamiltonian Eq. (|SJ) 
state given in Eq. ||SJ). Using the fact that the state in Eq. arises from successive Bogoliubov 
transformations one can calculate the expectation values. This is given as 

e = {n{(i,^J)\'HW,^J)) 

= /"dk|k|(2sin2/i^(k)+cos2/i,(k)(l-F"-Fi''^)) +yi + V2, (9) 

^ ' i,a 

where, the first term arises from the kinetic energy term and the second and third terms from the 
interation term of the Hamiltonian. Explicitly these are given as 

^i-yE fEc'-2E^r'), (10) 

i=l,2 \a=l,3 a=l,3 / 



with, 



= ^—^ / dk(i^*" - Fi") (11) 



(2^) 

related to vector density as p = (■0*"^' = 2 J^i a ■ Similarly 

" (2^ / '^^^^ " " ^i""^ ^'"^ '^^'^^^ ^^^^ 

is related to the scalar density — 2 a ■ '^^'^ effect of diquark condensates and/or density 
dependences are encoded in the functions i^*°(k) and F^'°(k) given as 

F*''(k) = sin^ 6'!_"(k) + sin^ f"" (e^^e"'' cos^ 6'i''(k) - sin^ 6l!?(k)) (l - 5""^) , (13) 

and, 

Fl^ild, ^ sin^ 6';°(k) + sin^ (cos^ 6'!^"(k) - e'^e'^'' sin^ fi*- (k)) (l - 5"^) . (14) 
Similarly the term V2 is given by 

V2 = -\g^llll (15) 

which arises from the diquark condensate terms only with 

II = j dk (sin2/*(cos2 9'} - sin^ 0^^\e,j\) + (cos^ 0'_l - sin^ 9^^\e,j\)^ cos{hi - hj) (16) 
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To calculate the thermodynamic potential we shall have to specify the chemical potentials 
relevant for the system. Here we shall be interested in the form of quark matter that might 
be present in compact stars older than few minutes so that chemical equilibriation under weak 
interaction is there. The relevant chemical potentials in this case then are the baryon chem- 
ical potential /is = 3fi, the chemical potential fiE associated with electromagnetic charge 
Q = (img(2/3, — 1/3, — 1/3) in flavor space, and the two color electrostatic chemical potentials 
/i3 and /ig corresponding to U{1)3 x U{l)s subgroup of the color gauge symmetry generated by 
cartan subalgebra Qs — diag(l/2, — 1/2, 0) and Qs — diag{l/3, 1/3, —2/3) in the color space. Thus 
the chemical potential is a diagonal matrix in color and flavor space, and is given by 

^'■^j,ab = {fJ-Sij + QijfJ'E)^ab + {Qsab + QsabfJ's)^ij ■ (17) 

Here, i,j are flavor indices and a, b are color indices. 
The thermodynamic potential is then given by, 



n = T + V- {fiN)-^s (18) 

where, we have denoted 

ifiN) = (^^'^Vy-.a^V-''") = 2 (19) 

i,a 

with being the chemical potential for the quark of flavor i and color a, which can be expressed 
in terms of the chemical potentials fi, /ZE,U3and fig, using Eq.^TJ. 
Finally, for the entropy density we have |l6| 

i.a 

+ sin^ e'_^ In sin^ +cos^d'_^ In cos^e'^y (20) 

Now if we minimise the thermodynamic potential with respect to hi (k) , we get 

, /, N (Mi — mi)k 
tan2/., k = . 21 

q + r Ea=l,3 Il '^i 



where. Mi — mi+ if' and = ^ (k^ + nif). Substituting this back in Ea. ll2|l we have the 

mass gap equation for the light quarks [nii—Q) 



(22) 



Clearly, the above includes the effect of diquark condensates as well as temperature and density 
through the functions F and Fx given in Eq.s H13|) and H14|l respectively. 

Next, minimising the thermodynamic potential with respect to the diquark condensate functions 
leads, to 

tan 2/2(k) = _ cos(/ii - hi) (23) 
e - V21 



tan 2/1 (k) = _ cos(/ii - /la) (24) 
e — j/ii 

In the above we have defined — (2(ji^/3)/||je,j |; e = (ei + 62) /2, Du = {vn + i/22)/2, V12 = 
(i^i2 + i^2i)/2. v^"" is the interacting chemical potential given as 

^za ^ _ (25) 

which may be expected in presence of a vectorial current-current interaction . In the above, 
p' ~ '^'Ylia=i 3 -^11° J with as defined in equantion Thus, it may be noted that the diquark 

condensate functions depend upon the average energy and the average chemical potential of the pair 
of quarks that condense. We also note here that the diquark condensate functions depends upon the 
masses of the two quarks which condense through the function cos(ft,i — h2)) with sin2/ii = Mi/ci, 
for u,d quarks which could be different when charge neutrality condition is imposed. Such a 
normalisation factor is always there when the condensing fermions have different masses as has 
been noted in Ref. Q| in the context of CFL phase. 

In an identical manner the di-antiquark condensate functions are calculated to be 

tan 2fl (k) = _ ^\ cos(/ii - /la) (26) 
e + vix 

tan 2/i2(k) = _ cos(/ii - /la) (27) 
e + V21 

Substituing the the solutions for the condensate functions in the definition for /| in Eq. H16(l . 
we have the superconducting gap equations given as 



£,"^21 + COs2(/li - h2) 
A2 



icos^ e^_} ~ Six? 



C^21 + Ai C0S2(/H - h2) 



1 /i21 • 2 /il2\ 

COS d, — sm d. 



cos2(/ii - /12) (28) 



and, 



A2 ^ r dk 

3(27r)3 ^ 



Ai 

f^^^+Af C0s2(/li - /l2) 



fcos2 0"-sin2 022) 



Ai 



^"2,, +A2C0S2(/11-/12) 



(cos^e^i-sin^ 022) 



,cos2(/ii - /i2) (29) 



In the above, t,±ia = e ± v^°' ■ Finally, minimisation of the thermodynamic potential with respect 
to the thermal functions ^?±(k) gives 

1 



sin^ 



exp(/3w5?) + 1 



(30) 



Various modes w*°s are given as follows. 



u?-} = w_ + ,5, - 8]} 
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UJ^3 = OJ^ 



J"} 



22 r 

Ll1_ = Ul^ — Of 



Lui^ =UJ+ + S^ + Si 



LU^J =LO++Sf+ Si 



Lol^ =UJ,-Sf- S}^ 



< =cj+-S,~ St' (31) 

and, finally, for the noncondensing colors cjj? — e'±i/*'^. As already mentioned, the first index refers 
to flavor and the second index refers to color. Here uj± — ^ IS? cos2(/ii — /12) + — (ei — e2)/2 

is half the energy difference of the two quarks which condense and e.g. S^} = {vn — i'22)/2, is 
half the difference of the chemical potentials of the two quarks which condense. Note that in the 
absence of imposing the charge neutrality condition all the four quasi particles will have the same 
energy uj^. Thus it is possible to have zero modes when charge neutrality condition is imposed 
depending both upon the values of S^ and S^,- So, although we shall have nonzero order parameter 
A, there will be fermionic zero modes or the gapless superconducting phase 

We shall 

discuss more about it section 4. 

Next, let us focus our attention for the specific case of superconducting phase and the chemical 
potential associated with it. First let us note that the diquark condensate functions depend upon 
the average of the chemical potentials of the quarks that condense. Since this is independent 
of /i3 we can choose /ia to be zero. In that case, = fi + 2/3/iB + /ig/x/S = /i^^ and also 
/i^^ — fj, — 1/3/iB + fis/VS = /^^^- This also means that the average chemical potential of both 
the condensing quarks are the same and is equal to p, = fi + l/6fiE + fJ-s/VS- Further S^^^ — 
fJ-Ef^ — (5^/4)(pi — /32)/2 = S^^^ = S^. In such a situation, Eq.(|2H|) and Ea.<|^ are both identical 
and hence we shall have only one superconducting gap equation given as 
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A = 4^ I dk 

3(27r)- 



— (COS^ „ gij^2 021\+A (cos2 9]^ ~ siu^ 9f] 



COS ( ^ ) (32) 



With this condition, it is clear from Eq.slPU that the quasi particle energies for each flavor becomes 
degenerate for both the colors which take part in condensation. Thus the quasi particle energies 
now become lo\ = lu- + S and UJ2 ~ — S, for u and d quarks respectively, with S ~ S^ — S^. 

Now making the use of this dispersion relations and the mass gap equations Ea. H22|) and the 
superconducting gap equation Ea. H32|l . the thermodynamic potential at zero temperature becomes. 
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with S — Se — S^, 



(27r)3 
4 



|k|--(c^_+^+) 



(27r)3 

3A2 "2 



5" 



(27r)3 



m - !le{k ^ kf) - 'jeim ~ kf) 



(33) 



The first three Hnes in Eq.JJJJ correspond to the contribution from the quarks taking part in the 
condensation while the last line is the contribution from the third color for the two light quarks. 

In an identical manner one can calculate the thermodynamic potential for the strange quark 
sector with strange quarks and anti-quarks for the vacuum structure. The contribution to the 
thermodynamic potential then is given by 

2 



^ ' 0=1,3 



|fc|2 + to2 _ 



k2 + Af326l(|k| - kf) 



(34) 



Here, pi — X]a=i 3 ^^"^ ~ 2/^,° for i, a = 1, 2, with given in Ea. (|ll|l . For the third color, 
— (^^43-j3^3^2 g^j^j jtqj, -j-j^g s^j-ange quark, p^a — (fc3'^)3/37r2. The fermi momenta are given by 
the usual relation k^^ = {v1^ - Mfy/"^. 

The total thermodynamic potential is given by including the contribution from the electrons and 
is given by 



127r2 



(35) 



Thus the thermodynamic potential is a function of four parameters: the three mass gaps and a 
superconducting gap which needs to be minimised subjected to the conditions of electrical and 
color neutrality. The electric and charge neutrality constraints are given respectively as 



^211 

Qe ^ gPi - gP2 - gP3 - Pe = 0, 



and, 



V3 



9*2 - 2p'^) = 0. 



(36) 



(37) 



Ea. (|35M57jl and the superconducting gap equation Ea. (|32|) constitute the basis of the numerical 
calculations that we discuss below. 

4. RESULTS AND DISCUSSIONS 



For numerical calculations we have taken the values of the parameters of NJL model as follows: 



g^A — 17.6, A — 0.68 GeV as typical values giving reasonable vacuum properties 



. We 
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FIG. 1: Masses of the quarks as a function of baryon chemical potential. Solid curve refers to masses of 
u and d quarks while the dashed curve refers to strange quark mass. 

might note here that the couphng introduced here is related to the usual scalar coupling of NJL 
model e.g. in Ref.|3 G as 5^ = SG. 

Current quark masses for u and d quarks are taken as zero and the current quark mass for strange 
quark is taken as 0.12 GeV. With this choice of parameters, the constituent quark masses at zero 
temperature and density are given as Mi = 0.35 GeV=Af2, and for strange quark M3 — 0.575 
GeV. 

In Fig^ we have plotted the variation of masses with baryon density without diquark conden- 
sation and without imposition of the charge neutrality conditions. The decrease of masses with 
density is a reflection of the decrease of quark condensates with density. 

Let us next discuss the case with diquark condensates along with quark antiquark condensates 
without imposing the charge neutrality condition. The numerical calculation for this case proceeds 
as follows. For a given chemical baryon potential, (and with ij^—O^ /is), the thermodynamic 
potential Ea. (|35|l is minimised with respect to the quark masses subjected to self consistently 
determining the interacting chemical potential using Eq. H25|) and solving the superconducting gap 
equation Ea. (|32|l . In Fig. 2 we have plotted the masses and the superconducting gap as a function 
of baryon chemical potential. The behaviour of masses are almost same to that without the 
diquark condensates. Here, however, we observe that for a small window in the quark chemical 
potential ( about 23 MeV), we have simultaneous existence of chiral symmetry breaking and 
color superconducting phase with the gap reaching a maximum upto 65 MeV in this window. 
Including vector interactions similar conclusion of simultaneous existence of both the condenstes 
was noted earlier in Ref. 9]. However, in Ref. the vector interaction did not contribute to the 
superconducting gap equation nor to the mass gap equation unlike the case here. We have compared 
the pressure in both the cases in this regime and it appears that in this small window, existence of 
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both the condensates has higher pressure than having only the quark antiquark condensates in the 
ground state. Beyond this window chiral symmetry is restored and the gap increases monotonically 
with chemical potential till the effect of cut off is felt and then it decreases. 
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FIG. 2: Masses of the quarks and the superconducting gap as functions of baryon chemical potential for 
the case fiE =0—)j,s- Solid, dashed and dot-dashed curves refer respectively to masses of up (down) quarks, 
strange quark and to the superconducting gap. 

We next discuss the results when charge neutrality conditions are imposed. As earlier, let us 
focus our attention first to the case without the diquark condensates. The results are shown in Fig. 
3. Here, the d- quark mass vanishes earlier than that of u quarks as the baryon chemical potential 
is increased. The reason is that to maintain electrical charge neutrality conditions the d quark den- 
sities are almost twice that of u quarks which makes the d quark antiquark condensates to vanish. 
In contrast to the rather sharp fall of the masses as compared to without imposition of electrical 
neutrality conditions, the u quarks remain massive much after d quark become massless (about 
80 MeV in the quark chemical potential window). The magnitude of electric chemical potential 
increases with densities till strange quarks begin to appear beyond which it starts decreasing so 
that charge neutrality conditions are maintained. 

Next we show the results with diquark condensates with neutrality conditions. For a given quark 
chemical potential /i, the interacting chemical potentials are determined using Ea. H25|) with a trial 
value of fiE and fi^- For high baryon chemical potential when chiral symmetry is restored for light 
quarks, the thermodynamic potential is varied with respect to the strange quark constituent mass 
after solving the superconducting gap equation Ea. (|32|l . The values of and /xg are varied so that 
the charge neutrality conditions Ea. H3t)|l and Ea. (|37|l are satisfied. The results are shown in Fig. 4. 
The behaviour of the superconducting gap is similar to that when charge neutrality condition are 
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FIG. 3: Masses of the quarks without superconducting phase when electrical charge neutrality condition is 
imposed. Solid, dot-dashed and dashed curves refer respectively to masses of up, down and strange quark 
mass. 




FIG. 4: Superconducting gap when color and electrical charge neutrality condition is imposed (solid curve). 
The dashed curve corresponds to when this condition is not imposed (i.e.fj,E=0). 

not imposed except that the magnitude of the gap decreases. The pressure of color and electrical 
neutral superconducting phase is lower than when neutrality conditions are not imposed as shown 
in Fig|Sl Further, the pressure of color and electric neutral normal quark matter is always lower 
than that of the color and electrically neutral superconducting matter. 

At high densities when chiral symmetry is restored for the light quarks, Se = and S = —Si, = 
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FIG. 5: Pressure as a function of baryon chemical potential when color and electrical charge neutrality 
condition is imposed (solid curve). The dashed curve corresponds to when this condition is not imposed. 
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FIG. 6: Number densities of u quarks (solid) and d quarks (dot-dashed) participating in superconducting 
phase. The density of s- quarks is also plotted (dashed curve). 



— (/k_e/2 — (5^/8)(p" — p'^)). For 5^ < 0, uji, the quasi particle energy for u quarks is always 
positive and hence only the d quark contribution will be there in the second term of the thermody- 
namic potential given in Ea. (|33|l . Further, it is easy to show that such a contribution arises when 
the magnitude of S^, which is half the difference of the chemical potentials of the quarks which 
condense, is the same or greater than the superconducting gap, A. When \6^\ — A, the mode 



LO2 — v^A^ + (e2 _ jy2^ _ §^ becomes gapless at the fermi sphere. For |i5y| > A, the gapless modes 
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occur at momenta higher than the (average) fermi momenta. In the present calculation this is the 
case for baryon chemical potential below 1600 MeV. The occurrence of such gapless superconduct- 
ing modes in neutral quark matter was first emphasised in Ref . [22! . Because of this the number 
densities of u and d quarks participating in the superconducting phase are not the same in this 
region. This is plotted in Fig. It may be noted that even in the gapless superconducting mode, 
the density of strange quarks is nonzero and can be as large as around 30-40 % of either of the 
condensing quarks. 

5. SUMMARY 

We have analysed here in a current current point interaction model, the structure of vacuum in 
terms of quark antiquark as well as diquark pairs. The methodology uses an explicit variational 
construct of the trial state and is not based on a mean field calculations. Because of the point 
interaction structure we could solve for the gap functions explicitly. The distribution functions 
are also determined variationally. Because of the vector interaction, the chemical potentials are 
interaction dependent and are calculated self consistently. We find that there is a small window in 
baryon chemical potential (about 80 MeV) when both the condensates are nonzero. 

To consider neutron star matter, we have imposed the condition of color and electric charge 
neutrality conditions through introduction of appropriate chemical potentials. It has been noted 
and emphasized earlier that projecting out the color singlet state from color neutral state costs 
negligible free energy for large enough chunk of color neutral matter JLOj . 

The gap reduces when color and electric neutrality conditions are imposed. The pressure with the 
gap is always higher than free quarks when charge neutrality conditions are imposed. For slightly 
lower densities, but large enough to be in the chiral symmetry restored phase, there appears to 
be gapless modes available when the condensing quarks have a difference in chemical potentials 
which is larger than twice the superconducting gap. In all these calculations we have also included 
the effect of self consistently determined mass for the strange quark. In fact, in the gapless 
superconducting phase, the number densities of strange quarks is also nonzero. 

We have focussed our attention here to the superconducting phase. The variational method 
adopted can be directly generalised to include color flavor locked phase and one can then make a 
free energy comparison regarding possibility of which phase would be thermodynamicaly favourable 
at what density. This will be particularly interesting for cooling of neutron stars with a CFL core. 
We have considered here homogeneous phase of matter. However, we can also consider mixed 
phases of matter with matter being neutral on the average. Some of these problems are being 
investigated and will be reported elsewhere [2^ . 
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